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Abstract 

■ Energy resolution, a/ (3 ratio, pulse-shape discrimination for 7 rays and a particles, 
radioactive contamination were studied with neodymium doped yttrium-aluminum garnet 
(YAG:Nd). Applicability of YAG:Nd scintillators to search for 2/3 decay and a activity of 

O ■ natural neodymium isotopes are discussed. 
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^ : 1 INTRODUCTION 

o : 

The great interest to the double beta (2/3) decay search is related, in particular, with the recent 
^ ' evidence of neutrino oscillations which strongly suggests that neutrinos have nonzero mass 
'q ■ [1, 2, 3, 4]. While oscillation experiments are sensitive to the neutrinos mass difference, only 

the measured neutrinoless (Ou) double beta decay rate could set the Majorana nature of the 

neutrinos and give the absolute scale of the effective neutrino mass. 
^ ■ 150 Nd is one of the most promising candidate for 2(3 decay study due to its large transition 

energy ((32/3) • As a result, the calculated value of the phase space integral G^ m of the 0u2(3 

■ decay of 150 Nd is the largest one among 35 possible 2(3~ decay candidates [5, 6]. The theoretical 
predictions of the product TfY 2 ■ (m u ) 2 are in the interval 3.4 x 10 22 - 3.4 x 10 24 yr (half-lives for 
the neutrino mass mechanism with (m u ) = 1 eV) [7]. Moreover, the larger the Q2/3 energy, the 
simpler, from an experimental point of view, is to overcome background problems. Note that 
the background from natural radioactivity drops sharply above 2615 keV, which is the energy 
of the 7 from 208 T1 decay ( 232 Th family). In addition, a contribution of cosmogenic activation, 
which is important problem of the next generation 2/3 decay experiments [8], decreases at higher 
energies. 

There exist no appropriate detector containing neodymium, which could serve as both source 
and detector simultaneously (so called "active" source experimental method [2]). It should be 
noted that this experimental approach will assure high registration efficiency which is especially 
important if very expensive enriched 150 Nd isotope would be used in an experiment. 

The purpose of our work was investigation of scintillation properties and preliminary check 
of radioactive contamination of neodymium doped yttrium-aluminum garnet (YAG:Nd) as a 
possible detector for double beta decay experiment with 150 Nd. The active source method allows 
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Table 1: Properties of YAG:Nd crystal scintillators 



Density (g/cm 3 ) 


4.56 


Melting point 


1970°C 


Crystal structure 


Cubic Garnet 


Hardness (Mohs) 


8.5 


Refractive index 


1.82 


Average decay time* 


4 /is 


Photoelectron yield relatively to Nal(Tl)* 


8% 



*For 7 rays, at room temperature. 



also to investigate other rare processes in Nd nuclei, such as a decay of naturally occurring 
neodymium isotopes. 

2 MEASUREMENTS AND RESULTS 

2.1 Energy resolution and a/j3 ratio 

Main properties of YAG:Nd (chemical formula YaAlsO^Nd) are presented in Table 1. These 
crystals are well known and their production is well developed due to their wide applica- 
tion in laser technique. The scintillation properties of cerium doped yttrium-aluminum garnet 
(YAG:Ce) have been studied in [9], properties of YAGYb have been reported in [10], whereas 
YAG:Nd, to our knowledge, was never studied as a scintillator. 

All the measurements were carried out for the 017 x 6 mm YAG:Nd crystal with ^2 
mol% of Nd. The mass of the crystal is 7.16 g. Photoelectron yield was estimated with 
the Philips XP2412 bialkali photomultiplier (PMT) as 8% of Nal(Tl). The energy resolution 
FWHM=13.6% was measured for 662 keV 7 line of 137 Cs with the YAG:Nd crystal wrapped by 
PTFE reflector tape and optically coupled to the PMT XP2412. A substantial improvement 
of the light output (~22%) and energy resolution was achieved by placing the crystal in liquid 
(silicone oil with index of refraction sal. 5). The crystal was fixed in center of the teflon container 
070 x 90 mm and viewed by two PMTs XP2412. Fig. 1 demonstrates the energy spectra of 
137 Cs and 207 Bi obtained in such a way with the YAG:Nd crystal. The energy resolution of 
9.3% was obtained with the YAG:Nd crystal scintillator for 662 keV 7 line of 137 Cs. 

We estimate the energy resolution of the YAG:Nd scintillation detector at the energy Q20 
of 150 Nd using the results of measurements with 137 Cs and 207 Bi 7 sources. These data were 
fitted by function FWHM(keV) =a + Jb x E y with values a = 2 keV and b = 5.2 keV, where 
energy of 7 quanta E y is in keV. The energy resolution FWHM^4% could be achieved at the 
Q 2 /3 energy of 150 Nd. 

The a/ (3 ratio was measured with the help of the collimated a particles of a 241 Am source. 
As it was checked by surface-barrier detector, the energy of a particles was reduced to about 
5.25 MeV by ~ 1 mm of air, due to passing through the collimator. Fig. 2 shows the energy 
spectrum of the a particles measured by the YAG:Nd scintillator. The a/ (3 ratio is 0.33. The 
YAG:Nd crystal was irradiated in three perpendicular directions with aim to check a possible 
dependence of the a signal on direction of irradiation. While earlier such a dependence has 
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Figure 1: Energy spectra of 137 Cs (a) and 207 Bi (b) 7 rays measured with the YAG:Nd scintil- 
lation crystal. The crystal was located in liquid and viewed by two distant PMTs (see text). 



been found for CdW0 4 [11] and ZnW0 4 [12] scintillators, we did not observe it for YAG:Nd 
crystal. 

2.2 Pulse-shape discrimination for 7 quanta and a particles 

The shapes of scintillation light pulses in YAG:Nd crystals were studied for 5.25 MeV a particles 
and ~ 1.5 — 1.8 MeV 7 quanta with the help of a 12 bit 20 MHz transient digitizer as described 
in [11, 13]. The pulse shape can be fitted by sum of exponential functions: 

/(*) = £ A/^ - r ) • ( e -*/n _ e -t/r 0): t > o, 

where are intensities (in %), and Tj are decay constants for different light emission compo- 
nents, t is integration constant of electronics («0.2 /is). The values of Ai and Tj obtained by 
fitting the average of ~3 thousand individual a and 7 pulses in the time interval — 60 /is (see 
Fig. 3) are presented in Table 2. Observed difference in light pulse shapes allows to discriminate 
7(/3) events from a particles. We applied for this purpose the optimal filter method proposed in 
[14] and developed in [13]. To obtain the numerical characteristic of YAG:Nd signal, so called 
shape indicator (SI), the following formula was applied for each pulse: 

SI = Ef(t k )-P(t k )/Ef(t k ), 
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Figure 2: The energy spectrum of a particles from collimated 241 Am source. 



where the sum is over time channels k, starting from the origin of pulse and up to 50 /xs, /(£&) 
is the digitized amplitude (at the time t k ) of a given signal. The weight function P(t) is defined 
as: P(t) = {f a (t) — fy(t)}/{fa(t) + f-y(t)}i where f a (t) and f 7 (t) are the reference pulse shapes 
for a particles and 7 quanta. Reasonable discrimination between a particles and 7 rays was 
achieved using this approach as one can see in Fig. 4 where the shape indicator distributions 
measured by the YAG:Nd scintillation crystal with a particles (E a xs 5.3 MeV) and 7 quanta 
(f« 1.5 — 1.8 MeV) are shown. 

2.2.1 Radioactive contamination of the YAG:Nd crystal 

To determine activity of a active nuclides from U/Th contamination in the YAG:Nd crystal, 
the pulse shape analysis was applied to data of 14.8 h low background measurement carried out 



Table 2: Decay time of YAG:Nd scintillator for 7 quanta and a particles. The decay constants 
and intensities (in percentage of the total intensity) are denoted as and Aj, respectively. 



Type of irradiation 


Decay constants, /is 


n (AO 


r 2 (A 2 ) 


r 3 (A 3 ) 


7 rays 

a particles 


1.1 (10%) 
1.0 (15%) 


4.1 (87%) 
3.9 (83%) 


20 (3%) 
18 (2%) 



4 




20 40 60 

Time (jis) 

Figure 3: Decay of scintillation in YAG:Nd for 7 rays and a particles (~ 3000 forms of 7 and a 
events were added) and their fit by sum of three exponential components with decay constants 
Tj « 1, « 4, and ss 20 /xs. 



at the earth surface laboratory in Kiev. The YAG:Nd crystal was viewed by the low-radioactive 
PMT (EMI D724KFLB) through plastic scintillator light-guide 10 cm in diameter and 4 cm 
long. The active light-guide reduces effect of 7 radiation from the PMT and provides suppression 
of cosmic rays induced background due to the pulse-shape discrimination of scintillation signals. 
Passive shield was consisted of steel (2 cm) and lead (10 cm). The energy resolution of the 
detector was 11% for 662 keV 7 rays. 

To estimate contamination of the crystal by a active nuclides of U/Th families, the a events 
were selected from the data accumulated with the YAG:Nd crystal with the help of the pulse- 
shape discrimination technique. In the a spectrum there is no peculiarities which can be surely 
attributed to a active U/Th daughters. Analyzing these data we set the limit on total a activity 
of 235 U, 238 U, and 232 Th daughters in the YAG:Nd crystal: < 20 mBq/kg. 

3 DISCUSSION 

3.1 2(3 decay of neodymium 

Among neodymium isotopes there are three possible 2(3 candidates: 146 Nd, 148 Nd, and 150 Nd. 
They are listed in Table 3 where the energies and the isotopic abundances are given. 
As it was mentioned in the Introduction, the most interesting of them is 150 Nd due to the 
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Figure 4: The shape indicator (see text) distributions measured by YAG:Nd scintillation crystal 
with a particles (E a = 5.25 MeV) and 7 quanta (m 1.5 — 1.8 MeV). 



large value of the transition energy. Two neutrino 2(3 decay of 150 Nd with the half-life Ti/ 2 = 
1-9-aI x 10 19 y r was observed in experiment [17] by using a time projection chamber and 
samples of enriched 150 Nd and natural neodymium. In [18] the double beta decay of 150 Nd 
was studied with similar technique (only enriched 150 Nd sample was used). The two-neutrino 
half-life T 1/2 = (6.8 ± 0.8) x 10 18 yr was measured. Recently, the observation of 150 Nd 2v2(3 
decay to the second excited level of 150 Sm (Of at 741 keV) with T 1/2 = 1.4lg;| x 10 20 yr was 
reported in [19], while only the limit Ti/ 2 > 1.5 x 10 20 yr was determined previously at 90% 
C.L. [20]. The most stringent limit on the neutrinoless mode 7i/ 2 > 1.2 x 10 21 yr was set in 
the experiment [18]. 

YAG:Nd scintillators provide a possibility to search for 0u2(3 decay of 150 Nd by "active" 
source experimental method. The studied pulse-shape discrimination ability is an important 



Table 3: 2(3 unstable neodymium isotopes 



Transition 


Isotopic abundance (%) [15] 


Q 2/3 (keV) [16] 


146 Nd^ 146 Sm 

148 Nd ^148g m 
150 Nd _^150g m 


17.2(0.3) 

5.7(0.1) 

5.6(0.2) 


70.2(2.9) 
1928.8(1.9) 
3367.5(2.2) 
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Table 4: Theoretical calculations of half-lives for a decay of natural Nd isotopes. Uncertainties 
given for calculations with cluster model of Ref. [25] and semiempirical formula of Ref. [26] are 
related only with uncertainty in the Q a values. For 143 Nd, we take into account that change in 
a parity and spin additionally suppresses the 143 Nd decay rate. 



Isotope Abund. Q a , MeV 
[15] [16] 


Calculated Ti/2, yr 
based on [25] based on [26] [27] 


Exp. Ti/2, yr 


143 Nd 12.2% 0.520(7) 
144 Nd 23.8% 1.905(2) 
145 Nd 8.3% 1.578(2) 
146 Nd 17.2% 1.182(2) 
148 Nd 5.7% 0.599(3) 


> l.ll 7 ; 8 xl0 80 5.2^ 7 xl0 93 
1.9 + [J- 2 xl0 15 4.3+^xlO 15 
1.7 + ° :2 xl0 22 3.9 + [J :3 xl0 23 3.7xl0 22 
2.0 + [J- 4 xl0 34 3.9 + °- 8 xl0 34 
6.1 + ^xl0 70 l.l+j^xlO 71 


2.29±0.16xl0 15 
> lxlO 17 



advantage of YAG:Nd scintillators as a low counting rate detector. Background of scintillation 
detector caused by radioactive contamination by U/Th daughters can be effectively rejected 
due to pulse-shape analysis as it was demonstrated in experiments with CdW04 [11, 21] and 
CaW0 4 [8, 22] crystal scintillators. Obvious disadvantage of YAG:Nd detector is rather small 
concentration of Nd. In the studied crystal the mass concentration is only ^0.03% of 150 Nd. It 
should be noted, however, that YAG:Nd crystals can be produced with neodymium concentra- 
tion up to 8 mol%, which corresponds to «0.1 mass% of 150 Nd. In this connection we would 
like to refer the project CANDLES [23] which intends to use a few tons of non-enriched CaF2 
crystals to search for 0z/2/3 decay of 48 Ca. The concentration of target 48 Ca nuclei in CaF2 
detector is at the same level: psO.1 mass% of 48 Ca. 

The value of energy resolution FWHM=4% could allow to discover [24] the 0z/2/3 decay of 
150 Nd with the half-life ~ 10 25 Y r ! "which corresponds to the effective Majorana neutrino 

mass m v ~ 0.06 — 0.6 eV (taking into account all existing calculations of nuclear matrix elements 
for 0u2(3 decay of 150 Nd). As for the sensitivity [24], an experiment involving tons of non- 
enriched YAG:Nd crystals with ps;8 mol% of Nd could reach the half-life sensitivity « 3 x 10 26 
yr (supposing zero background during ten years of measurements), which corresponds to the 
Majorana neutrino mass 0.01 — 0.1 eV. 

3.2 a decay of neodymium 

Alpha decay is allowed energetically for five naturally occurring isotopes of neodymium. Their 
energies of a decay (Q a ) and natural abundances are listed in Table 4. Calculation of half-life 
values based on models [25, 26], and the value of the half-life of 145 Nd calculated in [27] are also 
presented. Alpha decay was observed only for 144 Nd. Average half-life Ti/ 2 = 2.29 ±0.16 x 10 15 
yr [28] was derived on the basis of four measurements: T 1/2 = 2.2 x 10 15 yr [29], T 1/2 = 1.9 x 10 15 
yr [30], T 1/2 = 2.4 ± 0.3 x 10 15 yr [31] and T 1/2 = 2.65 ± 0.37 x 10 15 yr [32]. It is obvious that, 
due to low Q a values and long expected Ti/ 2 's, there is no perspective to observe the a decay 
of 143 Nd, 146 Nd and 148 Nd. As for 145 Nd, only two T 1/2 limits are known: T 1/2 > 6 x 10 16 yr 
[30] and T 1/2 > 1 x 10 17 yr [33]. 
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The indication on the a decay of natural 180 W isotope with half-life T±/ 2 = 1-2 jZoil (stat) 
±0.3 (syst) x 10 18 yr has been recently observed in the experiment [11] with the help of the low 
background 116 CdW0 4 crystal scintillators. This result was confirmed in the measurements with 
CaWC>4 crystal as scintillator [22] and bolometer [34]. It should be also referred the excellent 
result on the detection of alpha decay of 209 Bi with the half-life T ± / 2 = (1.9 ± 0.2) x 10 19 y 
[35]. In both the experiments [34] and [35] the cryogenic technique, which uses simultaneous 
registration of heat and light signals, have been applied. This method provides perfect selection 
of a events on the background caused by 7 rays (electrons). It should be interesting to check 
YAG:Nd crystal as the cryogenic detector to search for alpha decay of 145 Nd and to measure 
more accurately the 144 Nd half-life. 

It should be also mentioned a possibility to search for spin-dependent inelastic scattering 
of weakly interacting massive particles (WIMP) with excitation of low energy nuclear levels 
of 145 Nd (the lowest one: 3/2~ 67 keV, E2 transition, and the second level: 5/2" 72 keV, Ml 
transition) due to nonzero spin of this nucleus (7/2~). Identification of such "mixed" (nuclear 
recoil plus 7 quanta) events could be possible due to the simultaneous registration of heat and 
light signals (such a techniques have been used in [36, 37]). A heat /light ratio for such events 
would differ from "pure" nuclear recoils or 7(electron) events. 

4 CONCLUSIONS 

The scintillation properties of YAG:Nd crystals were studied. The energy resolution 9.3% 
(662 keV 137 Cs 7 line) was obtained with the YAG:Nd crystal scintillator placed in liquid and 
viewed by two PMTs. Shapes of scintillation signals were investigated, and reasonable pulse- 
shape discrimination for 7 rays and a particles was achieved. The a/ (3 ratio was measured with 
the YAG:Nd scintillator to be equal 0.33 at the energy of alpha particles 5.25 MeV. Radioactive 
contamination of the YAG:Nd crystal by a active nuclides from U/Th chains was estimated as 
< 20 mBq/kg. 

Three potentially 2(3 active neodymium isotopes can be studied with the help of YAG:Nd 
crystal. Due to good scintillation characteristics and pulse-shape discrimination ability the 
YAG:Nd scintillators seems to be encouraging material for the 2(3 experiment with 150 Nd. 

The YAG:Nd crystals could be used to search for a decay of natural neodymium isotopes, in 
particular, 145 Nd. A strong signature of a decay can be assured via simultaneous registration 
of heat and light signals using cryogenic technique. However, ability of YAG:Nd crystal as 
cryogenic detector has to be elaborated in an experiment. 
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